Mutant number distribution in an exponentially growing population

We present an explicit solution to a classic model of cell-population growth introduced by Luria and Delbrück [1] 70 years ago to study the emergence of mutations in bacterial populations. In this model a wild-type population is assumed to grow exponentially in a deterministic fashion. Proportional to the wild-type population size, mutants arrive randomly and initiate new sub-populations of mutants that grow stochastically according to a supercritical birth and death process. We give an exact expression for the generating function of the total number of mutants at a given wild-type population size. We present a simple expression for the probability of finding no mutants, and a recursion formula for the probability of finding a given number of mutants. In the “large population-small mutation”-limit we recover recent results of Kessler and Levin [2] for a fully stochastic version of the process.